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It  is  the  aim  of  the  present  work  to  address  some  of  the  aspects  of  microcracking  in  electron  beam 
deposited  scandia-stabilised  zirconia  electrolyte  applied  for  solid  oxide  fuel  cells  (SOFC)  where  a  thin 
electrolyte  layer  is  deposited  on  a  relatively  thick  anode  substrate.  A  model  of  microcracking  for  the  elec¬ 
trolyte  material  is  proposed  which  takes  into  account  the  statistical  distribution  of  grain  sizes,  the  stress 
redistribution  due  to  failure  of  individual  structural  elements  as  well  as  the  local  criterion  of  grain  frac¬ 
ture.  The  combination  of  electron  microscopy  research  with  model  calculations  permits  both  the  specific 
energy  of  new  surface  creation  in  the  electrolyte  and  critical  parameters  of  the  microcracking  process  to 
be  determined.  The  annealing-induced  electrolyte  microcracking  discussed  in  this  work  corresponds  to 
localised  microcracking,  where  each  next  structural  element  fails  mainly  at  an  existing  microcrack  tip. 
The  features  of  localised  microcracking  in  electron  beam  deposited  scandia-stabilised  zirconia  electrolyte 
are  analysed. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFC)  offer  promising  options  for  future 
energy  supply  due  to  the  high  electrical  energy  conversion  effi¬ 
ciency  they  display  (up  to  60%).  SOFC  have  a  particularly  wide 
range  of  applications  ranging  from  centralised  MW  scale  generation 
through  to  localised  applications  at  100+  kW  level  for  distributed 
generation  to  local  domestic  generation  on  the  1-101<W  scale  [1] 
and  as  Auxiliary  Power  Units  on  board  of  vehicles  (automotive, 
maritime  and  aerospace).  There  is  also  a  wide  range  of  other 
applications  where  cleaner  energy  is  required,  such  as  corrosion 
protection,  uninterruptible  power  supply,  remote  generation,  and 
domestic  appliances.  Another  recent  fuel  cell  development  is  the 
miniaturised  “palm-power”  units  of  lOOmW-lOW  for  portable 
electronic  devices.  Commercialisation  of  these  applications  will 
require  even  lower  operating  temperatures  and  additional  material 
improvements  [2]. 

A  single  SOFC  consists  of  an  air  electrode  (cathode),  electrolyte, 
and  a  fuel  electrode  (anode)  [3  ].  The  electrolyte’s  function  is  to  allow 
the  transport  of  oxide  ions  from  the  electrolyte’s  interface  with  the 
air  electrode  to  its  interface  with  the  fuel  electrode.  The  function  of 
the  air  electrode  is  to  facilitate  the  reduction  of  oxygen  molecules 
to  oxide  ions  transporting  electrons  to  the  electrode/electrolyte 
interface  and  allowing  gas  diffusion  to  and  from  this  interface.  The 
function  of  the  fuel  electrode  is  to  facilitate  the  oxidation  of  the  fuel 
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and  the  transport  of  electrons  from  the  electrolyte  to  the  fuel  elec¬ 
trode  interface.  In  addition  it  must  allow  diffusion  of  fuel  gas  to  and 
exhaust  gases  away  from  this  interface. 

Still,  several  major  development  problems  remain:  degradation 
is  too  high  (between  0.5  and  2%  per  1000  h  of  operation,  resulting  in 
a  ‘lifetime’  of  around  10,000-30,000  operating  hours)  [4];  costs  are 
too  high  (above  10,000€/kW);  intermittent  operation  (thermal  and 
oxidation  cycling)  results  in  increased  degradation  and/or  failure; 
power  density  at  temperatures  below  1023  K  is  comparatively  low 
[5].  These  problems  are  a  consequence  of  the  high  temperatures, 
elevated  pressures  required  for  their  operation,  and  the  ceramics 
materials  properties  [6]. 

The  use  of  electrolytes  with  higher  ion  conductivity,  for  instance 
scandia-stabilised  zirconia  (ScSZ),  offers  the  potential  for  lowering 
the  operating  temperature  with  an  ensuing  reduction  of  corrosion 
and  degradation  effects.  SOFC  operating  between  873  and  1073 1< 
are  generally  labelled  “intermediate  temperature”  SOFC  (IT-SOFC). 
On  the  other  hand,  the  higher  conductivity  can  also  be  read  as  a 
potential  for  obtaining  very  much  higher  power  densities  at  oper¬ 
ating  temperatures  above  1023  K.  I.e.,  these  materials  do  not  only 
offer  an  improvement  in  lifetime  of  SOFC,  they  could  also  be  used 
to  increase  the  compactness  of  Auxiliary  Power  Units  on  board  of 
vehicles  under  volume  and  weight  constraints.  Another  way  of  min¬ 
imising  ohmic  losses  across  the  electrolyte  is  to  use  thin  electrolyte 
membranes  with  thicknesses  between  5  and  50  pirn.  Electron  beam 
deposition  (EBD)  is  a  promising  technique  to  produce  thin  films  of 
electrolyte  on  thick  anode  substrates  of  planar  geometry.  The  anode 
supported  planar  SOFC  cell  has  become  one  of  the  most  widespread 
SOFC  types. 
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A  challenge  with  the  planar  geometries  is  in  obtaining  mechan¬ 
ically  stable  structures,  as  thin  layer  ceramics  are  inherently 
susceptible  to  failure  when  subjected  to  moderate  stress.  If  the 
anode  is  produced  from  a  nickel  oxide-ScSZ  composite  and  the  elec¬ 
trolyte  from  scandia-stabilised  zirconia,  thermal  mismatch  stresses 
keep  the  electrolyte  layer  under  compression  at  room  tempera¬ 
ture  and  under  tension  at  temperatures  higher  than  the  deposition 
temperature. 

The  mechanical  properties  of  the  SOFC  components,  including 
the  electrolyte,  must  be  known  in  order  to  design  an  optimal  SOFC 
structure.  Typically,  the  electrolyte  shows  microcracking  due  to 
residual  stress  [7-9].  Electrolytes  deposited  by  electron  beam  tech¬ 
nique  are  thin  films  with  specific  columnar  structure  attached  to 
a  weak  anode  layer.  It  is  impossible,  though,  to  obtain  bulk  speci¬ 
men  for  mechanical  testing  with  a  microstructure  which  is  identical 
to  the  electron  beam  deposited  electrolyte.  Furthermore,  inden¬ 
tation  methods  are  unable  to  determine  the  in-plane  electrolyte 
mechanical  properties.  Firstly,  unknown  residual  stresses  due  to 
edge  effects  change  the  indentation  response  [10];  secondly,  the 
columnar  microstructure  and  average  grain  size  comparable  with 
the  film  thickness  do  not  allow  imprints  which  contain  enough 
grains  (at  least  10)  for  the  necessary  averaging  of  properties.  In  such 
a  way,  the  determination  of  the  mechanical  properties  of  electron 
beam  deposited  electrolytes  is  complicated  scientific  problem. 

One  possible  approach  to  solve  the  problem  of  description 
and  prediction  of  the  mechanical  behaviour  of  such  objects  is 
microcracking  modelling.  Generalising  existing  structural  models 
of  mechanical  behaviour  of  materials  with  microcracks  [11-18],  it 
should  be  noted,  that  all  of  them  contain  the  following  issues:  struc¬ 
tural  element  (SE)  identification;  local  failure  criterion  for  SE;  stress 
redistribution  due  to  material  damage;  total  failure  criterion. 

It  is  the  aim  of  the  present  work  to  address  some  of  the  aspects 
of  microcracking  in  electron  beam  deposited  scandia-stabilised  zir¬ 
conia  electrolytes  applied  on  a  fuel  cell  anode.  In  order  to  gain 
better  insight  into  the  mechanical  behaviour  of  the  electrolyte,  the 
structural  model  of  scattered  microcracking  in  thin  brittle  films 
was  developed.  The  work  focuses  on  the  determination  of  both  the 
specific  energy  of  new  surface  creation  within  the  electrolyte  and 
critical  parameters  of  microcracking  processes  using  the  combina¬ 
tion  of  electron  microscopy  research  of  annealed  and  as-deposited 
electrolyte  with  model  calculations. 


slow  but  results  in  a  very  precise  control  of  layer  thickness  and, 
when  applied  expertly,  delivers  very  thin,  dense  coatings,  as  needed 
for  intermediate  temperature  electrolytes  [21].  These  thin  layers 
again  reduce  the  overall  resistance  of  the  electrolyte  and  further 
support  the  trend  towards  higher  power  density/lower  operating 
temperatures.  Powder  metallurgy  routes  were  applied  to  produce 
ScSZ  specimens  to  be  used  for  further  electron  beam  deposition 
experiments.  The  dimensions  of  the  cylindrically  shaped  specimens 
were  about  4  mm  in  diameter  and  about  8  mm  in  height  (porosity 
approx.  15%).  The  substrates  were  subjected  to  temperatures  of  873, 
1023,  and  1173  K  (deposition  temperature)  during  the  deposition 
process.  The  deposition  time  of  15  min  resulted  in  the  formation  of 
electrolyte  films  with  10  \xm  in  thickness. 

To  obtain  a  mechanical  loading  of  the  electrolyte,  half-cells 
were  annealed  at  1173,  1178,  1473,  1548  and  1633  K  for  2h.  Heat¬ 
ing  the  electrolyte  to  a  temperature  higher  than  the  deposition 
temperature  induces  equibiaxial  tension  due  to  a  thermal  expan¬ 
sion  mismatch  between  electrolyte  and  anode  layers.  The  thermal 
stresses  in  half-cells  during  annealing  can  be  calculated  from  layer 
properties  and  the  difference  between  annealing  and  deposition 
temperatures  using  the  method  presented  in  [22,23].  The  layer 
properties  required  for  the  model  calculations  are  elastic  moduli, 
Poisson  ratios,  coefficients  of  thermal  expansion  (CTE),  and  layer 
thicknesses. 

Standard  techniques,  scanning  electron  microscopy  and  Image- 
Pro  Plus  Software  (Version  3.0.00.00,  “Media  Cybernetics”),  were 
applied  to  obtain  electrolyte  microstructure  images  and  to  eval¬ 
uate  the  statistical  distribution  of  grain  and  microcrack  sizes  as 
well  as  the  number  of  microcracks,  average  microcrack  size  and  the 
standard  deviation  in  the  electrolyte  layer.  Using  Image-Pro  Plus 
Software  for  all  considered  states  of  electrolyte,  the  length  of  each 
microcrack  observed  in  corresponding  image  was  measured.  The 
microcrack  density  was  calculated  as 


(1) 


where  q  is  the  effective  length  of  the  ith  microcrack,  and  NA  is 
the  number  of  microcracks  on  the  scanned  area  A  (in  our  case 
A  =  45,000  |jim2).  The  number  of  microcracks  per  one  grain  n*/n* 
can  be  determined  from 


2.  Experimental 


n^  =  NA 

n|  Ng  ’ 


(2) 


The  composition  10mol.%  Sc203-1  mol.%  Ce02-89mol.%  Zr02 
(lOSclCeSZ)  was  found  to  be  optimal  for  the  electrolyte  and  as 
a  constituent  of  the  anode  [19,20].  Powders  manufactured  at  the 
Institute  for  Problems  of  Materials  Science  were  used  for  half-cell 
production  due  to  the  greater  sintering  ability  and  electric  conduc¬ 
tivity  compared  with  powders  from  other  manufacturers. 

SOFC  anodes  were  fabricated  with  a  traditional  ceramic  route 
[19].  Mixtures  of  40wt.%  NiO  and  60wt.%  lOSclCeSZ  were  ball- 
milled  with  alcohol  for  24  h.  Polyvinyl  alcohol  was  used  as 
pore-creating  agent.  Specimens  were  made  by  uniaxial  pressing 
using  pressure  ~20MPa.  These  were  sintered  at  1723 1<  for  2h  in 
air.  The  anode  substrates  were  found  to  meet  all  necessary  require¬ 
ments  concerning  both  strength  and  porosity  (100-150  MPa  and 
25-30%,  respectively).  The  thickness  of  the  anodes  was  990  p,m. 

Electron  beam  deposition  (EBD)  technique  was  used  for  produc¬ 
ing  electrolyte  films  deposited  on  the  porous,  non-reduced  anode 
substrates.  Hereafter  the  electrolyte-anode  layered  system  will  be 
termed  “half-cell”. 

EBD  is  one  form  of  physical  vapor  deposition  (PVD).  An  elec¬ 
tron  beam  is  directed  at  a  target  pellet,  in  order  to  ionize  and  melt 
the  target  material,  which  is  then  accelerated  in  an  electrical  field 
and  deposits  on  the  specimen  to  be  coated.  The  process  is  rather 


where  n*  is  the  number  of  cracks  per  area  unit,  n*  is  the  number  of 
grains  per  area  unit,  and  Ng  is  the  number  of  grains  on  the  scanned 
area  A. 

Thermal  expansion  analysis  of  anode  material  specimens  (of 
size  4  mm  x  4  mm  x  20  mm)  was  performed  to  obtain  information 
about  the  coefficients  of  thermal  expansion.  The  data  were  mea¬ 
sured  with  a  highly  sensitive  quartz  dilatometer  having  the  optical 
system  for  specimen  length  measurement.  Due  to  optical  magni¬ 
fication  within  the  range  4000 x  to  10,000 x  the  usage  of  such  a 
technique  allows  the  measurement  accuracy  of  specimen  length  to 
reach  0.00025-0.0001  mm.  CTE  measurements  were  made  in  air  in 
the  wide  temperature  range  of  293-1223  K.  The  heating  and  cool¬ 
ing  rates  were  10  and  20  K  min-1,  respectively.  The  quartz  thermal 
expansion  was  accounted  for  in  the  SOFC  constituent  CTE  measure¬ 
ments.  The  CTE  was  calculated  from  the  expression 


Q [t\-t2  — 


r(A2  -  A\ ) 
Rl(t2  -  ti) 


+  Oiq , 


(3) 


where  /  is  the  specimen  length;  r  is  the  roller  radius;  A\  and  A2  are 
the  specimen  length  values  corresponding  to  temperatures  t\  and 
t2,  respectively;  ft  =  1000  mm  is  the  roller-scale  distance;  aq  is  the 
CTE  of  quartz. 
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The  elastic  modulus  of  the  anode  material  was  determined  by 
using  four-point  bending  tests  with  an  Instron  testing  machine 
(internal  and  external  spans  were  20  and  40  mm,  respectively). 
Specimens  of  4  mm  x  5  mm  x  45  mm  in  size  were  used  for  the  tests. 

Disc-shaped  samples  of  anode  substrates  and  half-cells  were 
loaded  in  biaxial  ring-on-ring  bending  tests.  The  test  conditions 
followed  the  procedures  recommended  in  the  ASTM  standard  C 
1499-05  [24].  The  data  were  evaluated  using  the  equations  for  lin¬ 
ear  bending  theory.  Accordingly,  the  equibiaxial  flexure  stress  was 
derived  from: 


a  = 


3  P 
2ict2 


(4) 


where  P  is  the  applied  force,  t  is  the  specimen  thickness,  v  the  Pois¬ 
son  ratio  and  r\ ,  r2  and  r3  are  the  radii  of  the  load  ring,  supporting 
ring  and  (circular)  specimen,  respectively.  The  values  of  t,  ri,  r2 
and  r3  were  1,  9,  18.6  and  23.8  mm,  respectively.  Since  the  num¬ 
ber  of  tests  was  limited  (five  tests  for  anode  substrates  and  five 
tests  for  half-cells),  only  the  mean  value  and  the  standard  deviation 
are  given.  Equibiaxial  flexure  tests  were  carried  out  up  to  failure  at 
room  temperature  (RT)  in  air.  The  loading  rate  was  100  N  min-1  in 
all  tests. 


3.  Model 


The  structural  element  of  the  electrolyte  in  the  model  is  the 
individual  grain,  intact  or  cracked.  The  characteristics  of  the  indi¬ 
vidual  structural  elements  have  some  statistical  scatter  and  the 
loading  parameters  determine  the  number  of  cracked  elements. 
Each  structural  element  is  considered  as  located  in  effective  media. 
The  minimal  size  of  a  cracked  grain  and  the  failure  probability  can 
be  found  using  local  failure  criteria.  The  failure  probability  of  an 
individual  structural  element  will  determine  the  relative  number 
of  failed  elements. 

For  a  single-phase  electrolyte  the  main  parameter  in  determin¬ 
ing  the  stress  redistribution  is  its  elastic  modulus  E.  The  main 
loading  parameter  in  the  model  is  the  strain  energy  density: 


where  E  =  E(e),  i.e.  the  elastic  modulus  is  a  function  of  strain  e.  Addi¬ 
tionally,  the  elastic  modulus  of  an  undamaged  electrolyte  E0  also 
affects  the  stress  redistribution.  The  strength  of  an  individual  ele¬ 
ment  depends  on  the  specific  energy  of  the  new  surface  creation, 
y.  This  energy  is  the  surface  energy  of  the  material  in  the  simplest 
case  of  uniaxial  tension  and  transcrystalline  fracture. 

The  statistical  distribution  of  the  structural  element  sizes  is 
determined  by  few  factors  depending  on  the  material  fabrication 
technology.  Typically,  a  logarithmic  normal  distribution  of  grain 
sizes  is  used  [11,12,15,16].  Its  density  can  be  expressed  as: 


W)  = 


lDV2n 


exp 


/  (lni-JW)2\ 


(6) 


where  Z  is  the  grain  size;  M  is  the  mathematical  expectation  of  In  Z; 
D2  is  the  dispersion  of  the  distribution. 

It  should  be  noted  that  real  specimen  contains  finite  numbers 
of  grains.  Therefore,  if  grain  sizes  obey  the  logarithmic  normal  dis¬ 
tribution,  grains  with  size  larger  than  a  certain  size  Zmax  cannot  be 
observed.  The  number  of  grains  larger  than  Zmax  can  be  determined 
as  Nsp,  where  Ns  is  the  total  number  of  grains  in  the  specimen,  p  is 
the  probability  of  existence  of  a  grain  with  size  larger  than  Zmax  in 
the  specimen: 

pOO 

p=  m  di.  (7) 

•J  Jmax 


From  a  physical  point  of  view,  the  size  Zmax  can  be  defined  as 
that  size  for  which  the  total  number  of  grains  with  size  larger  than 
Zmax  equals  1.  In  fact,  this  results  in  absence  of  such  grains  in  the 
specimen.  In  such  a  way,  Zmax  can  be  determined  from  the  following 
equation  for  a  given  logarithmic  normal  distribution: 


Ns  /(/)  d/=  1.  (8) 

J  Imax 

A  strain  energy  density  criterion  is  used  as  the  local  failure  cri¬ 
terion  in  the  model.  The  failure  criterion  is  based  on  an  energetic 
approximation  of  crack  nucleation  conditions  [25].  It  allows  to  solve 
the  problem  of  a  random  flaw  in  a  grain  postulating  the  existence 
of  a  single  critical  defect  with  a  size  proportional  to  the  grain  size. 
The  local  failure  criterion  based  on  the  strain  energy  density  can 
be  expressed  as  fiq  >  qc,  where  qc  is  the  critical  value  of  the  strain 
energy  density  which  is  a  function  of  grain  size.  The  coefficient  /3 
characterises  the  strain  energy  density  and  stress  redistribution 
between  undamaged  material  and  the  effective  continuum  con¬ 
sisting  of  undamaged  and  cracked  structural  elements.  According 
to  [26],  j3  is  considered  to  be  approximated  by  the  ratio  of  elas¬ 
tic  moduli  of  undamaged  and  damaged  materials,  i.e.  in  our  case 
P  =  E0IE. 

The  parameter  qc  for  brittle  material  can  be  determined  from 
energy  balance  condition  as  follows: 


(9) 


where  I<  is  the  geometrical  factor  of  the  structural  element.  In 
many  important  applications  of  microcracking  models  we  have  a 
so-called  2D  case  when  a  thin  layer  of  material  contains  microc¬ 
racks  through  the  full  thickness  [27,28].  More  specifically,  a  thin 
layer  can  be  defined  as  the  layer  with  thickness  less  than  one  or 
two  average  grain  sizes.  In  this  case,  for  an  in-plane  isotropic  layer 
in  first  approximation,  the  individual  grains  can  be  considered  as 
circular  cylinders  with  height  equals  to  the  layer  thickness.  Then, 
from  the  ratio  of  crack  surface  area  and  structural  element  volume 
the  geometrical  factor  equals  8 In  for  transcrystalline  cracks  when 
a  crack  is  located  along  the  diameter  of  a  cylinder.  This  factor  is  8 
for  intercrystalline  cracks  when  the  crack  coincides  with  the  lateral 
surface  of  a  cylinder. 

The  dependencies  of  elastic  modulus  and  Poisson  ratio  on  the 
damage  parameter  (hereinafter  referred  to  as  microcrack  density) 
for  the  2D  case  presented  in  [27]  are: 


(10) 

(ID 


where  ff=n*(c 2)  is  the  microcrack  density,  c  is  the  crack  length, 
angular  parentheses  denote  averaging  over  all  cracks,  u0  is  the  Pois¬ 
son  ratio  of  undamaged  material. 

The  minimal  size  of  a  failed  grain  Z*  can  be  determined  from  the 
equation 


Pq  =  qc 


(12) 


and  Eq.  (9)  for  1  =  1*.  In  the  framework  of  our  failure  criterion,  all 
grains  with  size  larger  than  Z*  will  be  cracked  and  grains  with  size 
less  than  Z*  remain  undamaged.  Note  that  Z*  depends  on  microcrack 
density  in  material  through  /3. 

The  normalised  strain  energy  density  is  defined  as 

<?*=<?  y,  (13) 

where  lm  is  the  average  grain  size.  The  normalised  stress  and  strain 
are  a*  =  cr^lm/E0y  and  e *  =  e^E0lm/y,  respectively. 
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In  the  2D  case  the  microcrack  density  of  the  damaged  layer  can 
be  determined  as: 


W 


c2/c(c)  Ac, 


(14) 


where  (c2)  =  J  c2/c(c)dc,  /C(c)  is  the  density  of  the  statistical 
distribution  of  microcrack  sizes.  Under  the  condition  that  each 
microcrack  occupies  only  a  single  grain  (so-called  scattered  micro¬ 
cracking)  and  the  microcrack  size  equals  the  grain  size,  the  density 
of  the  microcrack  size  distribution  for  grains  with  size  larger  than 
P  can  be  found  from  the  grain  size  distribution: 


MC)  =  ~Ft, 


f(c) 


J,,m  di' 


(15) 


The  denominator  in  (15)  supplies  the  normalisation  of  the  dis¬ 
tribution  density  function.  Note  that/c(c)  =  0  for  cracks  with  size 
less  than  P  and  grains  with  size  larger  than  /max  are  not  considered. 
The  fraction  of  failed  grains  (microcrack  number  per  grain)  is 

Wmax 

m  di.  (is) 


*  /**  r 

n*  _  / 

"py,* 


Note  that 

Imax 


T  =  t  /  PfW  dl 

n*  4j 
0 


(17) 


is  the  average  area  of  a  single  grain,  if  we  consider  grains  as  circular 
cylinders  for  simplification. 

The  main  model  equations  for  scattered  microcracking  of  a 
single-phase  thin  layer  of  brittle  material  taking  into  account  (9), 


(10).  (12),  (15),  (16)  are: 

4  Ji!max/2/(/)  dl 
tf  ~  it 

(18) 

p  2KyE0 

a2  exp(7r///2)  ’ 

(19) 

f,imax  (f(l)dl 

|c,~  Jlrmti’ 

(20) 

_2  flre/wi  2 

J)r7mi  ’ 

(21) 

where  a  is  the  stress  acting  in  the  electrolyte;  (c)  is  the  average 
microcrack  size;  and  s  is  the  standard  deviation  of  microcrack  size. 
Using  substitution  of  (19)  to  (18),  the  dependence  a  on  ff  can  be 
calculated.  Taking  into  account 


Q  = 


2  E  2En  6XP 


(22) 


the  dependence  of  strain  energy  density  on  microcrack  density  can 
be  determined  as  well. 


4.  Results  and  discussion 

The  electrolyte  layer  in  an  IT-SOFC  is  a  single-phase  anisotropic 
material  with  the  anisotropy  resulting  from  its  processing  route. 
The  thin  layer  of  electrolyte  has  a  columnar  structure  where  grains 
are  orientated  perpendicular  to  the  specimen  surface  (Fig.  1  a).  Nev¬ 
ertheless,  such  a  structure  is  isotropic  in  the  plane  of  the  electrolyte, 
i.e.  the  properties  of  the  film  in  different  directions  within  the 
plane  are  identical.  The  main  structural  elements  of  the  electrolyte 
which  are  subject  to  microcracking  are  its  grains.  The  electrolyte 
microcracking  is  in  fact  a  2D  case  because  all  cracks  formed  in  the 
electrolyte  layer  are  oriented  perpendicular  to  its  plane.  According 


Fig.  1.  Microstructure  of  as-deposited  EBD  lOSclCeSZ  electrolyte  film  (deposition 
temperature  1173  K):  (a)  side  view  (intercrystalline  fracture  and  columnar  structure 
of  electrolyte)  and  (b)  top  view. 


to  SEM  investigations  the  intercrystalline  fracture  mode  is  revealed 
(Fig.  la). 

The  microstructure  of  the  as-deposited  electrolyte  film  (deposi¬ 
tion  temperature  1173  K)  is  shown  in  Fig.  lb,  where  the  boundaries 
between  the  structural  elements  are  distinctly  visible.  It  was  found 
experimentally  that  the  distribution  of  grain  sizes  in  the  electrolyte 
demonstrates  a  slight  dependence  on  deposition  temperature. 
Therefore,  in  first  approximation,  it  is  reasonable  to  use  a  common 
distribution  which  is  the  same  for  all  deposition  temperatures  in  the 
model.  The  distribution  of  grain  sizes  in  the  electrolyte  is  shown  in 
Fig.  2a,  the  distribution  of  normalised  grain  sizes  is  shown  in  Fig.  2b. 
According  to  experimental  data,  the  average  grain  size  in  the  elec¬ 
trolyte  is  6.64  |xm,  the  standard  deviation  is  1.53  p,m.  It  was  found 
also,  that  the  experimental  statistic  distribution  of  normalised  grain 
sizes  can  be  approximated  with  sufficient  accuracy  by  a  logarith¬ 
mic  normal  distribution  (Fig.  2b)  with  parameters  M= -0.02601, 
D2  =  0.052028,  /max  =  17  jxm. 

One  of  the  possible  approaches  to  characterise  the  mechanical 
properties  of  EBD  electrolyte  films  is  to  use  heat  treatment,  i.e. 
annealing  at  various  temperatures  above  deposition  temperature 
followed  by  electron  microscopy  investigation  of  the  electrolyte 
surface.  The  annealing  creates  equibiaxial  stress  states  with  tensile 
stress  in  the  electrolyte,  whose  value  depends  on  the  difference 
AT  between  annealing  temperature  and  deposition  temperature. 
Electron  microscopy  investigations  of  the  electrolyte  surface  have 
allowed  the  density  of  microcracks  arising  due  to  tensile  residual 
stresses  in  the  annealed  electrolyte  to  be  determined  (Fig.  3). 

The  initial  data  for  thermal  stress  calculation  in  half-cells  were 
established  experimentally  and  via  literature  survey.  The  elastic 
modulus  for  the  anode  material  was  found  to  be  81  GPa  from 
four-point  bending  tests.  The  elastic  modulus  of  an  undamaged 
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Fig.  2.  Statistical  distribution  of  structural  element  sizes  in  as-deposited  EBD 
lOSclCeSZ  electrolyte  film:  (a)  non-normalised  and  (b)  normalised;  the  curves  indi¬ 
cate  the  corresponding  logarithmic  normal  distributions. 


electrolyte  was  adopted  to  be  200  GPa  [29].  It  is  known  that  the 
elastic  modulus  of  ceramics  as  a  rule  demonstrates  a  dependence 
on  temperature.  For  example,  the  existence  of  such  dependence 
for  yttria-stabilised  zirconia  (YSZ)  is  shown  in  [30]  where  a  sub¬ 
stantial  scattering  of  experimental  data  is  observed.  The  scattering 
complicates  the  determination  of  a  reliable  dependence  of  elastic 
modulus  on  temperature.  To  the  best  of  our  knowledge,  the  tem¬ 
perature  dependence  of  the  elastic  modulus  for  lOSclCeSZ  is  not 
available  from  literature.  Therefore,  as  a  first  approximation,  a  con¬ 
stant  value  of  the  elastic  modulus  is  used  for  calculation  in  this 
work.  The  Poisson  ratios  of  anode  and  electrolyte  were  adopted  to 
be  0.3  and  0.31,  respectively  [30].  The  coefficient  of  thermal  expan¬ 
sion  for  the  anode  was  found  experimentally  to  be  11.7  x  10-6  K-1. 
The  coefficient  of  thermal  expansion  for  the  electrolyte  was  adopted 
to  be  1  x  10-5 1<-1  [31  ].  Note  that  non-reduced  anode  substrate  was 
used  in  this  work.  In  such  a  way,  there  are  no  particles  of  pure  nickel 
in  the  anode.  This  means  that  creep  processes  in  the  anode  during 
annealing  were  negligible,  simplifying  residual  stress  calculation. 

The  validity  of  the  initial  data  and  thermal  stress  calculation  pro¬ 
cedures  in  half-cell  specimen  is  confirmed  by  results  of  ring-on-ring 
mechanical  tests.  The  fracture  stress  of  anode  specimens  is  found  to 
be  41.5  ±4.0  MPa.  The  fracture  stress  of  half-cell  specimens  corre¬ 
sponding  to  anode  failure  is  28.1  ±3.0  MPa.  The  resulting  decrease 
of  anode  strength  (13.4  MPa  on  the  average)  is  associated  with 
residual  tensile  stresses  in  the  anode  near  the  anode/electrolyte 
interface.  The  calculated  value  of  residual  stress  in  the  anode  is 
16  MPa  (tension)  which  is  in  good  accord  with  the  anode  strength 
decrease,  confirming  hereby  the  validity  of  the  residual  stress  cal¬ 
culation. 

It  is  found  that  the  microcrack  density  is  negligible  at  A T=  150  K 
demonstrating  the  practical  absence  of  microcracking  due  to 
stresses  associated  with  this  temperature  difference.  Increasing  the 
temperature  difference  gives  rise  to  a  growth  of  the  microcrack 
density  which  promotes  the  significant  decrease  of  the  electrolyte 
elastic  modulus  and  Poisson  ratio  [27].  In  turn,  this  results  in 


Fig.  3.  Microstructure  of  EBD  lOSclCeSZ  electrolyte  film  after  deposition  and  annealing:  (a)  deposition  at  1173  K,  annealing  at  1473 1<;  (b)  deposition  at  873  K,  annealing  at 
1178  K;  (c)  deposition  at  1023  K,  annealing  at  1473  K;  (d)  deposition  at  873  K,  annealing  at  1473 1<. 


M.  Lugovy  et  al.  /  Journal  of  Power  Sources  194  (2009)  950-960 


955 


Table  1 

The  effective  elastic  properties  of  lOSclCeSZ  electrolyte  film  and  thermal  stresses 
corresponding  to  different  AT. 


decreasing  effective  thermal  stresses  in  the  electrolyte.  The  effec¬ 
tive  elastic  properties  of  the  electrolyte  and  acting  thermal  stresses 
corresponding  to  different  AT  are  listed  in  Table  1 .  The  dependence 
of  the  tensile  stresses  on  microcrack  density  can  be  compared  with 
the  model  calculations. 

The  results  of  statistical  analysis  of  microcrack  length  (his¬ 
tograms  and  logarithmic  normal  distributions)  for  different  AT  of 
300,  305,  375,  450,  460  and  600  K  are  shown  in  Fig.  4.  The  aver¬ 
age  microcrack  size  and  the  standard  deviation  which  correspond 
to  the  temperature  difference  AT=300K  are  11.79  and  5.56  pm, 
respectively  (Fig.  4a);  the  number  of  microcracks  revealed  exper¬ 
imentally  in  the  scanned  area  is  104.  At  a  temperature  difference 
AT=305  K  the  average  microcrack  size  and  the  standard  deviation 
equal  12.6  and  7.54  pm,  respectively  (Fig.  4b);  the  corresponding 
number  of  microcracks  is  99  in  this  case.  The  average  microc¬ 
rack  size  and  the  standard  deviation  corresponding  to  temperature 
difference  AT=375K  are  20.40  and  10.03  pm,  respectively;  the 
number  of  microcracks  is  118.  The  average  microcrack  size  and 
the  standard  deviation  corresponding  to  temperature  difference 
AT=450  K  are  26.53  and  10.95  pm,  respectively  (Fig.  4c);  the  num¬ 
ber  of  microcracks  is  115.  At  a  temperature  difference  AT = 460  K  the 
average  microcrack  size  and  the  standard  deviation  equal  25.58  and 
13.88  |jim,  respectively;  the  number  of  microcracks  revealed  is  128. 
At  a  temperature  difference  AT = 600  K  the  average  microcrack  size 


Fig.  5.  Logarithmic  normal  distributions  of  grain  sizes  and  crack  sizes  in  the  elec¬ 
trolyte  film:  1— grain  sizes;  2— crack  sizes  at  temperature  difference  of  300  K; 
3— crack  sizes  at  temperature  difference  of  305  K;  4— crack  sizes  at  temperature 
difference  of  450  K;  5— crack  sizes  at  temperature  difference  of  600 1<. 

and  the  standard  deviation  equal  29.29  and  15.79  pm,  respectively 
(Fig.  4d);  the  number  of  microcracks  revealed  is  137.  The  experi¬ 
mental  data  obtained  indicate  the  fact  that  increasing  AT  brings 
about  an  insignificant  growth  of  the  number  of  microcracks  which 
in  fact  can  be  considered  as  a  constant  value  within  the  statisti¬ 
cal  scatter  of  data.  At  the  same  time,  a  marked  growth  of  both  the 
average  microcrack  size  and  the  standard  deviation  are  observed. 
For  comparison,  the  appropriate  logarithmic  normal  distributions 
of  microcrack  sizes  along  with  the  logarithmic  normal  distributions 
of  grain  sizes  are  shown  in  Fig.  5.  It  is  seen  that  the  distribution  for 
microcrack  sizes  are  wider  compared  to  that  for  grain  sizes.  In  addi¬ 
tion,  the  distribution  of  microcrack  sizes  becomes  wider  with  AT 
increasing. 


c,  pm  c,  pm 


Fig.  4.  Histograms  and  the  corresponding  logarithmic  normal  distributions  of  microcrack  sizes  in  the  electrolyte  film  after  deposition  and  annealing:  (a)  deposition  at  1173  K, 
annealing  at  1473  K;  (b)  deposition  at  873  K,  annealing  at  1178  K;  (c)  deposition  at  1023  K,  annealing  at  1473  K;  (d)  deposition  at  873  K,  annealing  at  1473  K. 
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It  is  important  to  note  that  the  combination  of  experimental 
observations  of  microcrack  density  at  different  annealing  temper¬ 
atures  with  model  calculations  allows  the  critical  parameters  of 
electrolyte  failure,  such  as  maximum  strain  energy  density  Qc  and 
strain  energy  density  corresponding  to  microcracking  initiation  q0, 
as  well  as  the  corresponding  stresses  to  be  determined. 

As  a  rule,  the  calculation  of  microcracking  parameters  using 
(18)  and  (19)  gives  the  dependencies  of  stress  and  strain  energy 
density  on  microcrack  density  having  maxima  and  the  non-zero 
strain  energy  density  corresponding  to  zero  microcrack  density. 
This  means  that  a  certain  value  of  strain  energy  density  q0  is 
necessary  for  initiating  microcracks.  The  normalised  value  of  this 
parameter  q*0  is  given  by  expression  (13),  the  corresponding  stress 
value  is: 

(70  =  \/2<Jo  E.  (23) 

The  values  of  maximum  stress  and  normalised  maximum  stress 
are  designated  as  crmax  and  <7j^ax,  respectively.  The  microcrack 
density  corresponding  to  the  maximum  stress  is  less  than  that  cor¬ 
responding  to  the  maximum  strain  energy  density.  After  reaching 
maximum  strain  energy  density  the  microcracking  process  devel¬ 
ops  catastrophically  when  fracture  becomes  localised  and  each  next 
grain  fails  at  the  tip  of  the  existing  crack.  In  such  a  situation  the  total 
number  of  cracks  remains  practically  constant.  Eventually,  the  total 
failure  of  a  material  occurs  in  this  stage  of  microcracking.  Therefore, 
the  maximum  value  of  strain  energy  density  Qc  (the  normalised 
value  of  this  parameter  is  designated  as  Q^  and  given  by  ( 13 ))  can  be 
used  as  a  critical  parameter  corresponding  to  failure.  Note  that  the 
proposed  model  which  describes  scattered  microcracking  allows 
the  start  point  of  localised  microcracking  to  be  determined. 

In  case  of  uniaxial  tension,  the  critical  stress  oc  which  corre¬ 
sponds  to  strain  energy  density  Qc  is  as  follows 

0c  =  y/2 qJe.  (24) 

Note  that  for  brittle  materials  the  uniaxial  tension  strength 
equals  the  equibiaxial  tension  strength  [32].  This  means  that  in 
the  equibiaxial  stress  state  the  critical  stress  oc  can  be  set  to  be 
comparable  with  the  thermal  residual  stress  ar  in  the  electrolyte. 

The  parameter  necessary  for  model  calculations  but  unknown 
at  this  stage  is  the  specific  energy  of  new  surface  creation  y.  An 
attempt  is  made  here  to  determine  this  parameter  for  the  elec¬ 
trolyte  material  analysing  the  experimental  data  from  the  point  of 
view  of  some  limitations  imposed  on  microcracking  parameters. 
First  of  all,  the  dependence  of  the  average  microcrack  size  (c)  on  the 
minimal  size  of  failed  element  P  can  be  calculated  using  (20),  i.e. 
based  on  the  assumption  of  scattered  microcracking.  The  results  of 
this  calculation  are  shown  in  Fig.  6a,  curve  1.  In  this  case  all  points 
of  the  curve  (c)  =f[l*)  have  no  dependence  on  y.  Next,  each  of  the 
experimentally  determined  average  microcrack  sizes  can  be  consid¬ 
ered  as  being  in  correspondence  with  a  given  parameter  P  which 
can  be  found  from  (19)  prescribing  some  specific  value  to  y.  The 
position  of  a  point  corresponding  to  maximum  strain  energy  den¬ 
sity  (i.e.  to  localised  microcracking  initiation)  depends  on  y  too, 
but  in  fact  this  dependence  is  negligible.  Note  that  all  points  of 
curve  1  to  the  left  of  the  indicated  point  (which  corresponds  to 
localised  microcracking  initiation),  are  shown  for  comparison  only 
and  do  not  describe  any  real  process.  Because  curve  1  demonstrates 
the  decrease  of  parameter  (c)  upon  decreasing  P,  but  experimen¬ 
tal  points  at  the  temperature  difference  >300  K  for  any  value  of  y 
demonstrate  opposite  tendency,  i.e.  growth  of  the  average  micro¬ 
crack  size  (c)  upon  P  decreasing,  the  conclusion  can  be  made  that 
annealing  heat  treatments  caused  the  localised  microcracking.  Tak¬ 
ing  into  account  that,  in  line  with  the  model  assumptions,  localised 
microcracking  initiation  corresponds  to  achieving  maximum  strain 
energy  density  Qc,  and  a  fracture  is  an  irreversible  process  run- 


Fig.  6.  Dependencies  of  the  average  crack  size  (a)  and  microcrack  density  (b)  on 
minimal  size  of  failed  grain  for  electrolyte  film  lOSclCeSZ:  1— scattered  microcrack¬ 
ing;  2— localised  microcracking;  the  arrow  indicates  the  transition  from  scattered 
to  localised  microcracking;  filled  circles,  open  circle,  triangle  and  rhomb  desig¬ 
nate  experimental  points,  maximum  strain  energy  density,  maximum  stress  and 
microcracking  initiation,  respectively. 

ning  with  continuous  decrease  of  P,  it  follows  that  all  experimental 
points  in  Fig.  6a  may  not  be  located  to  the  right  of  the  point 
which  corresponds  to  Qc.  First  of  all,  this  restriction  concerns  the 
experimental  point  with  minimal  average  size  of  microcrack.  Thus, 
it  follows  from  (19)  that,  other  conditions  being  equal,  y  cannot 
exceed  a  given  value  yu;  in  other  words,  we  have  a  boundary  for 
parameter  y:  y<yu- 

Shown  in  Fig.  6b  is  the  dependence  of  the  microcrack  density 
ff  on  the  parameter  P  obtained  using  (18).  It  is  seen  that  for  a 
given  ff  value  the  points  corresponding  to  localised  microcrack¬ 
ing  (curve  2)  are  on  the  right  of  the  points  which  correspond  to 
scattered  microcracking  (curve  1).  According  to  (19),  other  condi¬ 
tions  being  equal,  the  points  for  localised  microcracking  correspond 
to  less  stress  as  compared  to  points  for  scattered  microcracking. 
Thus,  the  experimental  points  for  localised  microcracking  must  be 
located  lower  than  curve  1  describing  the  dependence  of  stress  on 
microcrack  density  for  scattered  microcracking  (Fig.  7).  Note  that 
the  decrease  of  the  parameter  y  will  shift  points  of  curve  1  down¬ 
wards.  In  other  words,  parameter  y  cannot  be  less  than  some  given 
value  yh  signifying  the  lower  boundary,  y>y[.  Based  on  the  exis¬ 
tence  of  the  above  two  restrictions,  in  our  case  the  lower  boundary 
yi  and  upper  boundary  yu  of  parameter  y  were  found  to  equal  0.035 
and  0.037Jm-2,  respectively.  This  means  the  best  estimation  of 
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Fig.  7.  Dependence  of  stress  on  microcrack  density  for  a  lOSclCeSZ  electrolyte 
films:  1— scattered  microcracking:  2— localised  microcracking:  the  arrow  indicates 
the  transition  from  scattered  to  localised  microcracking:  filled  circles,  open  circle, 
triangle  and  rhomb  designate  experimental  points,  maximum  strain  energy  density, 
maximum  stress  and  microcracking  initiation,  respectively. 


the  specific  energy  of  new  surface  creation  (with  accuracy  of  3%) 
is  y  =  0.036 Jnrr2. 

The  question  arises  to  what  measure  the  obtained  value  of  the 
specific  energy  of  new  surface  creation  is  realistic  and  how  it  com¬ 
pares  with  literature  data.  According  to  the  authors’  knowledge, 
the  surface  energy  of  lOSclCeSZ  has  not  been  determined  to  date. 
Experimental  values  of  the  surface  energy  of  pure  Zr02 ,  according  to 
different  references,  are  as  follows:  1.294Jm-2  (high-temperature 
wetting  experiments  with  liquid  metals)  and  0.0559Jm-2  (low 
temperature  wetting  experiments  with  biological  liquids)  [33]; 
6.5  ±  0.2  J/m2  (monoclinic  zirconia),  2.1  ±  0.05  J  nrr2  (tetragonal  zir- 
conia),  0.5  ±  0.05 Jnrr2  (amorphous  zirconia)  [34]— measured  by 
oxide  melt  solution  calorimetry;  0.77  and  1.13  J  nrr2  (tetragonal  and 
monoclinic  zirconia)  [35].  For  polycrystalline  Zr02  containing  8% 
Y203,  Tsoga  and  Nikolopoulos  [36]  extrapolated  high-temperature 
data  to  get  a  surface  energy  of  1.927  J  nrr 2  at  0  K.  It  is  obvious  that  the 
value  of  the  specific  energy  of  new  surface  creation  (0.036 J  nrr2 )  for 
lOSclCeSZ  electrolyte  obtained  in  our  work  is  reduced  essentially 
in  comparison  with  specific  surface  energy  values  of  zirconia-based 
ceramics  given  in  literature.  The  main  reason  of  the  reduction  in  our 
opinion  can  be  attributed  to  the  intercrystalline  fracture  mode  of 
electrolyte  film  (Fig.  1).  It  is  known  that  the  energy  of  intergran¬ 
ular  fracture  of  polycrystalline  materials  decreases  to  some  extent 
compared  to  the  energy  of  transcrystalline  fracture.  The  decrease  is 
widely  considered  to  be  associated  with  the  occurrence  of  impurity 
atoms  and  their  segregation  to  grain  boundaries.  The  segregation 
is  known  to  affect  both  mechanical  properties  and  conductivity  of 
ceramics  [37-42].  In  addition  to  this  possible  impurity  effect,  in  our 
opinion,  a  second  reason  for  the  low  specific  energy  of  new  surface 
creation  of  electrolyte  film  cannot  be  ruled  out.  This  is  associated 
with  an  intrinsic  weakness  of  grain  boundaries  of  the  electrolyte 
film  resulting  from  EBD  processing  and  from  its  columnar  struc¬ 
ture.  Considering  above  arguments,  it  is  evident  that  the  energy  of 
free  surfaces  formed  upon  intercrystalline  brittle  failure  should  be 
used  in  our  model  calculations  instead  of  the  surface  energy. 

Taking  the  specified  value  of  the  specific  energy  of  new  surface 
creation  as  0.036J  nrr2,  the  dependence  of  the  average  microcrack 
size  on  the  minimal  size  of  failed  grain  /*  was  found  to  be  described 
satisfactorily  by  the  expression 

(C>  =30-0.029/*3,  (25) 

where  (c)  and  P  are  in  micrometers.  Expression  (25)  indicates  an 
existence  of  ultimate  average  microcrack  size  being  30  pan.  Mean¬ 


while,  we  take  into  account  the  fact  that  a  microcracking  process 
runs  from  large-sized  to  fine-sized  grains  to  finish  theoretically 
with  zero-sized  grains. 

Taking  the  specified  value  of  y,  elastic  modulus  of  undamaged 
electrolyte  of  200  GPa  and  experimental  grain  size  distribu¬ 
tion,  the  critical  parameters  of  electrolyte  microcracking  were 
calculated  using  expressions  (18)  and  (19).  The  strain  energy 
density  corresponding  to  microcracking  initiation  and  maxi¬ 
mum  strain  energy  density  were  established  to  be  q0  =  0.017  MPa 
and  Qc  =  0.027  MPa,  respectively.  The  stress  which  corresponds 
to  microcracking  initiation,  maximum  stress  and  stress  corre¬ 
sponding  to  localised  microcracking  initiation  were  found  to  be 
or0  =  82MPa,  ormax  =  100MPa  and  ac  =  95  MPa,  respectively.  The 
microcrack  density  corresponding  to  maximum  stress  and  to 
localised  microcracking  initiation  were  found  to  equal  /max  =  0.077 
and  fc  =  0.226,  respectively. 

Taking  into  account  the  above  critical  parameters  and  using  the 
dependencies  of  elastic  modulus  and  Poisson  ratios  on  microcrack 
density  (10),  (11 ),  the  values  for  the  temperature  difference  AT  cor¬ 
responding  to  microcracking  initiation,  the  maximum  stress  and 
the  maximum  strain  energy  density  were  found.  The  temperature 
differences  are  195,  257  and  283  K,  respectively.  The  temperature 
differences  used  in  this  work  were  150,  300,  305,  375,  450,  460, 
and  600  K.  Therefore,  it  is  evident  that  the  temperature  difference 
of  150  K  has  not  caused  any  microcracking,  while  temperature  dif¬ 
ferences  of  300,  305,  375,  450,  460,  and  600  K  were  more  than 
283  K,  thus  corresponding  to  localised  microcracking.  Note  also 
that  the  experimental  determination  of  microcrack  density  in  the 
electrolyte  in  conditions  being  within  the  range  from  microcrack 
initiation  to  maximum  stress  is  a  complicated  task  due  to  the  fact 
that  the  corresponding  microcrack  densities  are  not  far  from  the 
so-called  background  value  which  is  observed  in  practically  non- 
cracked  material.  Thus,  it  is  realistic  for  scattered  microcracking 
to  reveal  the  conditions  after  maximum  stress,  i.e.,  being  within 
the  range  from  maximum  stress  to  maximum  strain  energy  den¬ 
sity.  In  our  case  this  is  the  range  of  temperature  difference  from 
257  to  283  K  which  is  a  narrow  temperature  interval  of  26  K  only. 
Accounting  also  for  the  occurrence  of  statistical  deviations  in  the 
microstructure  of  the  electrolyte  film,  it  is  a  difficult  task  to  fit  the 
annealing  conditions  into  the  temperature  range  of  distinct  scat¬ 
tered  microcracking. 

The  determination  of  parameter  y  was  based  on  the  model 
assumption  that  localised  microcracking  begins  at  the  moment  of 
achieving  maximum  strain  energy  density.  Fig.  8  is  the  additional 
evidence  in  favor  of  this  assumption.  Fig.  8a  shows  the  dependence 
of  a  number  of  cracks  per  one  grain  n*/n*  upon  microcrack  den¬ 
sity.  It  is  seen  that  experimental  values  n*/n*  for  all  microcrack 
densities  are  located  near  the  line  depicting  the  calculated  num¬ 
ber  of  cracks  per  one  grain  which  correspond  to  the  maximum 
strain  energy  density.  Taking  into  account  the  statistical  character 
of  experimental  values  n*/n*,  it  is  realistic  to  consider  that  in  first 
approximation  the  number  of  cracks  at  the  stage  of  localised  micro¬ 
cracking  remains  constant.  The  validity  of  the  model  calculations  is 
confirmed  also  by  Fig.  8b  where  both  theoretical  (scattered  microc¬ 
racking)  and  experimental  (localised  microcracking)  dependences 
of  the  average  microcrack  size  on  microcrack  density  are  shown. 
It  is  seen  that  these  dependencies  seem  to  intersect  in  the  point 
which  correspond  to  maximum  strain  energy  density.  In  addition, 
the  conclusion  can  be  made  that  the  average  microcrack  size  for 
localised  microcracking  increases  upon  growth  of ff,  demonstrating 
the  tendency  for  saturation. 

The  standard  deviation  s  of  microcrack  sizes  increases  while 
the  average  microcrack  size  (c)  decreases  in  the  stage  of  scattered 
microcracking  (Fig.  9,  curve  1 ).  Development  of  the  localised  micro¬ 
cracking  causes  a  further  increase  of  standard  deviation  and  growth 
of  the  average  microcrack  size  (Fig.  9,  curve  2).  The  experimen- 
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Fig.  8.  Dependencies  of  the  number  of  microcracks  per  one  grain  (a)  and  aver¬ 
age  microcrack  size  (b)  on  microcrack  density  for  lOSclCeSZ  electrolyte  films: 
1— scattered  microcracking:  2— localised  microcracking:  the  arrow  indicates  the 
transition  from  scattered  to  localised  microcracking:  filled  circles,  open  circle,  tri¬ 
angle  and  rhomb  designate  experimental  points,  maximum  strain  energy  density, 
maximum  stress  and  microcracking  initiation,  respectively. 


Standard  Deviation,  pm 


Fig.  9.  Dependence  of  average  microcrack  size  on  its  standard  deviation  for 
lOSclCeSZ  electrolyte  films:  1— scattered  microcracking:  2— localised  microcrack¬ 
ing:  the  arrow  indicates  the  transition  from  scattered  to  localised  microcracking: 
filled  circles,  open  circle,  triangle  and  rhomb  designate  experimental  points, 
maximum  strain  energy  density,  maximum  stress  and  microcracking  initiation, 
respectively. 


Fig.  10.  Dependence  of  critical  microcracking  parameters  on  the  parameter  y/ys 
for  lOSclCeSZ  electrolyte  films:  (a)  strain  energy  densities,  (b)  critical  stresses; 
1— maximum  values;  2— values,  corresponding  to  microcracking  initiation;  I— zone 
of  undamaged  material  (below  curve  2);  II— zone  of  scattered  microcracking 
(between  curves  1  and  2);  III— zone  of  localised  microcracking  (above  curve  1);  the 
dotted  line  is  the  boundary  between  transcrystalline  and  intercrystalline  fracture. 


tal  dependence  (i.e.,  localised  microcracking)  of  (c)  upon  s  can  be 
expressed  in  first  approximation  as  (c)  =  2s. 

It  is  of  interest  to  know  in  what  measure  a  variation  of  the  param¬ 
eter  y  will  affect  the  critical  parameters  of  the  electrolyte  film.  The 
calculations  were  made  taking  the  value  of  the  true  surface  energy 
of  the  electrolyte  film  in  first  approximation  to  be  ys  =  1.927  J  m-2 
[36].  The  dependence  of  the  normalised  strain  energy  density  on 
parameter  y\ys  for  an  electrolyte  film  of  10  pm  in  thickness  is 
shown  in  Fig.  10.  Curve  1  is  the  maximum  strain  energy  density 
which  corresponds  to  failure;  curve  2  is  the  strain  energy  den¬ 
sity  which  corresponds  to  microcracking  initiation.  There  are  three 
characteristic  zones  of  the  material  behaviour.  Zone  I  (below  curve 
2)  is  the  zone  of  existence  of  undamaged  material;  zone  II  (between 
curves  1  and  2)  and  zone  III  (above  curve  1 )  corresponds  to  scattered 
microcracking  and  catastrophic  localised  microcracking,  respec¬ 
tively.  Note  that  on  the  left  of  the  dotted  line  the  material  will 
be  susceptible  to  intercrystalline  fracture.  In  this  region  the  strain 
energy  density  is  a  linearly  increasing  function  of  the  parameter 
ylys.  On  the  right  of  the  dotted  line  lies  the  range  of  transgran- 
ular  failure  where  the  strain  energy  density  does  not  depend  on 
ylys.  According  to  the  model  assumptions,  the  critical  value  of 
y\ys ,  corresponding  to  the  transition  from  intercrystalline  to  tran¬ 
scrystalline  fracture,  is  the  ratio  of  the  geometrical  factors  of  the 
structural  elements  for  transcrystalline  cracks  to  that  for  intercrys¬ 
talline  cracks,  which  in  the  2D  case  is  0.3183.  The  dependencies 
of  stress  corresponding  to  microcracking  initiation  (curve  2)  and 
maximum  stress  (curve  1)  on  y\ys  are  shown  in  Fig.  10b.  Similarly 
to  Fig.  10a,  three  different  zones  can  be  found.  The  dotted  line  is 
the  boundary  between  transcrystalline  and  intercrystalline  frac- 
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ture.  Within  the  zone  of  intercrystalline  fracture  the  critical  stresses 
increase  with  growth  of  the  parameter  y/ys,  but  these  dependen¬ 
cies  are  non-linear.  Within  the  zone  of  transcrystalline  fracture  the 
critical  stresses  do  not  depend  on  y\ys . 

The  second  issue  to  be  discussed  is  the  effect  of  the  electrolyte 
layer  thickness  on  the  microcracking  process.  It  is  obvious  that  the 
layer  thickness  will  affect  the  microcracking  density  corresponding 
to  saturation.  The  average  distance  between  microcracks  is  lim¬ 
ited  by  the  layer  thickness  in  having  an  effect  on  the  saturation 
microcracking  density.  In  our  case,  the  minimum  average  distance 
between  microcracks  (for  annealing  with  AT=600K)  is  not  less 
than  20  p>m  (a  factor  two  higher  than  the  layer  thickness).  Actu¬ 
ally,  we  only  observed  a  tendency  to  saturation  but  the  saturation 
condition  itself  was  not  reached.  In  [43]  the  calculation  procedure 
for  critical  stress  (depending  on  layer  thickness)  to  maintain  tunnel 
crack  growth  was  proposed.  It  was  found  that  this  is  conservative: 
any  flaw,  regardless  of  its  initial  size  or  shape,  cannot  grow  into  a 
long  tunnel  if  the  applied  stress  is  below  the  critical  stress.  Using 
the  parameters  of  the  lOSclCeSZ  electrolyte  the  critical  stress  for 
crack  tunneling  was  calculated  to  be  45  MPa  which  is  less  than  the 
critical  stresses  calculated  using  the  proposed  model.  This  means 
that  the  occurrence  of  microcracking  in  the  lOSclCeSZ  electrolyte 
is  also  in  accord  with  [43].  In  addition,  the  proposed  model  takes 
the  layer  thickness  indirectly  into  account.  According  to  the  model, 
all  critical  parameters  of  the  microcracking  process  depend  on  the 
average  grain  size.  However,  in  the  2D  case  considered  in  this  work, 
the  average  grain  size  is  comparable  with  the  layer  thickness. 

5.  Conclusions 

The  combination  of  electron  microscopy  data  with  model  cal¬ 
culations  permits  both  the  specific  energy  of  new  surface  creation 
of  the  electrolyte  and  critical  parameters  of  the  microcracking  pro¬ 
cess  to  be  determined.  The  best  estimation  of  the  specific  energy  of 
new  surface  creation  (with  an  accuracy  of  3%)  is  y  =  0.036 Jrrr2. 
The  strain  energy  density  of  microcracking  initiation  and  maxi¬ 
mal  strain  energy  density  are  0.017  and  0.027  MPa,  respectively. 
The  mechanical  stress  corresponding  to  microcracking  initiation  is 
82  MPa;  the  maximum  mechanical  stress  is  100  MPa;  the  mechan¬ 
ical  stress  corresponding  to  initiation  of  localised  microcracking 
is  95  MPa.  The  microcrack  densities  which  correspond  to  maximal 
stress  and  to  the  initiation  of  localised  microcracking  are  0.077  and 
0.226,  respectively. 

The  annealing-induced  electrolyte  microcracking  described 
here  corresponds  to  localised  microcracking,  when  each  next  struc¬ 
tural  element  fails  mainly  at  existing  microcrack  tips.  The  following 
features  of  localised  microcracking  in  EBM  deposited  lOSclCeSZ 
electrolyte  are  revealed:  (i)  the  average  microcrack  size  increases 
with  microcrack  density  growth,  demonstrating  a  tendency  for  sat¬ 
uration;  (ii)  the  standard  deviation  of  microcrack  size  is  about  half 
of  the  average  microcrack  size;  (iii)  the  number  of  microcracks 
increases  to  a  little  degree;  in  first  approximation  it  is  a  constant 
corresponding  to  a  microcrack  number  at  a  moment  when  strain 
energy  density  is  at  a  maximum.  It  is  found  also  that  the  min¬ 
imum  value  of  average  microcrack  size  is  reached  at  a  moment 
corresponding  to  maximum  strain  energy  density. 

However,  it  should  be  noted  that  the  problem  of  electrolyte 
microcracking  is  far  from  a  comprehensive  understanding.  Among 
the  most  prospective  directions  of  future  investigations  the  issues 
of  further  collection  of  experimental  data  as  well  as  development  of 
theoretical  foundations  of  localised  microcracking  can  be  identified. 
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